Zero-temperature equation of state of quasi-one dimensional H2 
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We have studied molecular hydrogen in a pure ID geometry and inside a narrow carbon nanotube 
by means of the diffusion Monte Carlo method. The one-dimensionality of H2 in the nanotube is 
well maintained in a large density range, this system being closer to an ideal ID fluid than liquid ^Ife 
in the same setup. H2 shares with ''He the existence of a stable liquid phase and a quasi-continuous 
liquid-solid transition at very high linear densities. 



PACS numbers:67.70.+n,02.70.Lq 

The experimental finding in 1991 [|l| of carbon nan- 
otubes opened brand new possibilities both in technol- 
ogy and in fundamental physics. The nanoscale of these 
new materials has led to the discovery of novel mechan- 
ical, chemical and electrical properties [|| that suggest 
exciting new technological applications. One of the more 
relevant features of these new materials is their large 
adsorption capability compared with that of a graphite 
planar substrate. In fact, the enhanced Van der Waals 
interaction when particles are adsorbed inside a single 
nanotube, or in the interstitial sites of single wall carbon 
nanotube (SWCN) bundles, may allow the density of the 
adsorbed substrate to be high enough for the bundles to 
serve as gas storage devices. Special interest exists in the 
physisorption of hydrogen (^-|7j in the quest for a fuel 
cell efficient enough to be used as a pollution-free energy 
carrier. It has been recently suggested that SWCN with 
diameters of the order of a nanometer can be the best 
candidates to approach the pursued level of packing . 

From a more fundamental point of view, the strong 
confinement of particles adsorbed in the carbon chan- 
nels of a SWCN bundle, with diameters ranging form 7 
to 40 A and an width to length ratio of ^ 1000, offers 
the possibility of an experimental realization of a quasi- 
one dimensional system. Moreover, if the temperature is 
low enough one is dealing with a unique opportunity of 
studying a nearly one-dimensional quantum fluid. In a 
recent experiment, Teizel et al. have unambiguously 
observed the quasi-one dimensional behavior of ^He ad- 
sorbed in a SWCN bundle by measuring its desorption 
rate. On the other hand, theoretical studies in the limit of 
zero temperature and strictly one dimension have proved 
the existence of a liquid state with a binding energy in 
the milli-Kelvin scale [p|-p]|. Using the diffusion Monte 
Carlo (DMC) method, we also compared the ideal ID 
geometry system with a realistic situation in which ^He 
is adsorbed in a narrow nanotube jl^. The aim of the 
present work is to extend the theoretical study to the 
appealing case of H2 that, besides its technological rele- 
vance, might offer the existence of a homogeneous liquid 
phase at zero temperature. It is worth noticing that both 



bulk and two-dimensional II2, are solid in this temper- 
ature limit. Liquid phases have only been observed in 
theoretical calculations of small clusters jlj] and in 2D 
geometries with localized alkaline impurities . 

We have studied molecular hydrogen at zero temper- 
ature in a one-dimensional (ID) array and inside a sin- 
gle walled carbon nanotube (T) of radius R = 3.42 A, 
which corresponds to a (5,5) armchair tube iQ. The 
technique used is the DMC method, which has become 
in the last decades one of the most efficient theoretical 
tools, from the microscopic point of view, to deal with 
quantum fluids. The DMC method solves the iV-body 
imaginary-time Schrodinger equation 



dt 
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in which /(R,t) i/'(R)4'(R,i), D = /2m, ElCR) = 
'0(R)^-^i/V(R) is the local energy, and F(R) — 
2 -(/^(R)"^ Vi/'(R) is the so-called drift force. The time- 
independent wave function ipCR) acts as an important 
sampling function and is chosen as a trial model usually 
adjusted at a variational Monte Carlo (VMC) level. In 
bosonic systems like the present one the DMC method 
provides exact results apart from statistical uncertain- 
ties. More specific details of DMC may be drawn from 
Refs. |5|jl|. 

A relevant issue in a microscopic study is the nature of 
the interspecies interaction. We have considered the II2 
molecules interacting via the isotropic semiempirical po- 
tential from Silvera and Goldman (SG) ||l^ that has been 
extensively used in path integral Monte Carlo (PIMC) 
and DMC calculations of bulk ||l8| , clusters ^ and II2 
films |l|. The SG is a pair potential that incorporates 
to some extent the effect of three-body interactions by 
means of an effective two-body term of the form Cg/r^. 
On the other hand, the isotropy of the potential is well 
justified if one considers that at very low temperatures al- 
most all the H2 molecules are para- hydrogen species, i.e., 
they are in the J=:0 rotational state. In the simulations 
of H2 inside a nanotube, we consider a cylindrically sym- 
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metric potential as suggested by Stan and Cole iQ. In 
that simplified model, the interactions between C atoms 
and H2 molecules are axially averaged out resulting in 
a potential which only depends on the distance to the 
center of the tube. It has been proved |^ that the differ- 
ences between that smoothed potential and a potential 
which is built up as an explicit sum of individual C-II2 
interactions are not significant and surely smaller than 
the relative uncertainty in the (<T,e) Lennard- Jones pa- 
rameters. Considering a — 2.97 A and e ~ 42.8 K, the 
symmetric potential felt by a II2 molecule in a (5,5) tube 
has a depth of 42 e, three times larger that the attraction 
of the same molecule in a flat graphitic surface. 

The use of importance sampling in DMC requires the 
introduction of a trial wave function that guides the dif- 
fusion process to relevant regions of the walkers phase 
space. In the ID system, we consider 



with «'j(R) = n»<j exp 



2 (rt ) 
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a Jastrow wave 



function with a McMillan two-body correlation factor. 
Inside the nanotube, II2 molecules interact with the walls 
of the cylinder and therefore we have added an additional 
one-body term 



(3) 



with 'I'c(R) = cxp{—crf), being the radial dis- 
tance of the particle to the center, to avoid the hard core 
of the H2-nanotube interaction. 

Theoretical calculations of ID "'He agree in predict- 
ing a liquid-solid phase transition at high linear densities 
[p|-pl|. This transition which is only possible at abso- 
lute zero temperature looks like a nearly continuous one 
without a measurable difference between the melting and 
freezing densities. Following the same procedure that in 
our previous work in helium, we have explored the exis- 
tence of such a transition in II2. In this ordered phase, 
we modify the trial wave function in both the ID sys- 
tem and inside the nanotube by multiplying them by a 
z-localized factor 5's(R) = J|f^ exp(— a(zi — ZigY). The 
sites Zis are equally-spaced points in both the ID line 
and the axial direction of the nanotube. 

The variational parameters a, b and c have been op- 
timized by means of VMC calculations. In the liquid 
phase and near the equilibrium density b— 3.759 A, with 
a slight increase with density (at A — 0.277 A^^, b ~ 
3.789 A), whereas c = 4.908 A'^ is kept fixed for aU A 
values. In the sohd phase, h = 3.404 A, a = 0.799 A"^ 
and c — 5.136 A~^, with a negligible A dependence in the 
region analyzed. 

The possible existence of a liquid-solid phase transi- 
tion at high linear density has been studied in both ID 
and inside a narrow nanotube. In Table I, results for 
the energy per particle of both systems are reported for 



the liquid (a = 0) and solid (a ^ 0) phases. The com- 
parison between the energies of both phases at the same 
density shows that their difference changes sign in going 
from A = 0.312 k'^ to A = 0.304 k'^ in ID and from 
A = 0.320 A-i to A = 0.312 A-^ in the tube. Above 
these densities, the system prefers to be localized in a 
solid-like structure with a difference \E{s) — E{1)\ that 
increases with A. When the density decreases the liquid 
phase is energetically preferred and again the size of the 
difference \E{s) — E{l) \ increases when A diminishes. The 
density value at which this difference becomes zero is es- 
timated to be A = 0.309 A^^ in ID and A = 0.315 A^^ in 
the tube, being not possible to distinguish between freez- 
ing and melting densities. As previously studied in ''He 
|lO| it appears to be a nearly continuous phase transition 
located at a density close to the inverse of the location 
of the minimum of the respective pair potential {r^ — 

0. 337 A-i vs \s = 0.358 k'^ for helium, and r^^ = 0.291 
A"-' vs As = 0.309, 0.315 A~^ for molecular hydrogen). 

Inside the nanotube (T), the energies are much more 
negative that in ID due to the strong attraction of the 
carbon substrate: the binding energy of a single H2 
molecule in the tube is E^ = —1539.87 ±0.11 K. Looking 
at the T-energy results contained in Table I one realizes 
that also in this case a transition occurs at a density 
very close to the ID one. It is remarkable that both in 
ID and T, H2 remains liquid below the liquid-solid transi- 
tion density, and thus a homogeneous liquid phase at zero 
pressure is predicted. That result contrasts with the the- 
oretically and experimentally well established solid phase 
in 3D and the 2D solid phase predicted by a PIMC 
calculation p9| ]. 

The equations of state of liquid H2 near the equilibrium 
density for both the ID and T systems are shown in Fig. 

1. In order to make the energy scales compatible we have 
subtracted the single binding energy Ef, to the T results. 
The lines in the figure correspond to the third-degree 
polynomial fits in the form 



B 



An 



(4) 



The best set of parameters Cq, Aq, A and B are re- 
ported in Table II. The equilibrium densities in both sys- 
tems are the same considering their respective uncertain- 
ties but the binding energy eg = e(Ao) is larger when H2 
is inside the nanotube. The difference between the 3D ge- 
ometry (T) and the idealized one (ID) can be quantified 
by means of the adimensional parameter 
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In the present system, around Aq, A"^ = 3.5 % which 
emphasizes the proximity between the real system and 
the idealized one. It is worth noting that in ^He inside 
the same nanotube we obtained A"^ = 90 % [p^, and 
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therefore H2 seems a better candidate to experimentally 
achieve a ID condensed phase. That significant difference 
between helium and hydrogen may be understood taking 
into account that the difference between the hard-core 
size of the C-He interaction {ac-He = 2.74 A) and the 
C-H2 one {ac-H2 = 2.97 A) is magnified in a (5,5) tube 
because its small radius {R = 3.42 A). 

In Fig. 2, the density dependence of the pressure for 
both the ID and T systems is reported from equilibrium 
up to the liquid-solid transition density. As a matter of 
comparison, the same results for ^He are also plotted. 
Both in H2 and *He the pressure increases faster in the 
ID geometry (P\) than in the tube (P) due the trans- 
verse degree of freedom that particles have in the latter 
case (notice the proportionality between the scales of P 
and Pa in Fig. 2, Px/P = ttR"^). At a given density A, 
the difference between the T and ID pressures is smaller 
in H2 than in *He. For example, at the respective tran- 
sition densities that difference is more than one and a 
half times larger in helium than in hydrogen. There- 
fore, the one-dimensionality of H2 inside the nanotube is 
well maintained in all the liquid regime in contrast with 
^He, in which the departure from such an idealized model 
already appears around the equilibrium density and in- 
creases significantly with A. Also apparent form Fig. 2 
is a much smaller compressibility in H2 than in *He. In 
the ID geometry at A = Aq the velocity of sound in H2 is 
c = 736.1 ± 0.2 m/sec to be compared with c = 7.98 ± 
0.07 m/sec in ^He at A = AoC'He) = 0.062 A-^. The ve- 
locity of sound drops to zero at the spinodal point that, 
according to the equation of state (Table II), is located 
at densities 0.210 ±0.001 A^^ and 0.209 ±0.001 A^^ for 
the ID and T systems, respectively. 

The spatial structure of molecular hydrogen in the ID 
array and inside the nanotube has also been analyzed 
by means of the two-body radial distribution function, 
g(z), and its Fourier transform, the static structure fac- 
tor S{k). In Fig. 3, ID results for g{z) arc reported at 
both the equilibrium density for the liquid phase, and in 
the solid-liquid transition region (A = 0.312 A^^). The 
corresponding results for H2 inside the nanotube arc in- 
distinguishable from the ID results in the scale; shown 
in Fig. 3. The solid g{z) shows a strongly localized or- 
der around the equally spaced z-sites that decreases very 
slowly when z increases. The result for g{z) at A Aq 
manifests the nature of a dense fluid with an appreciable 
structure that decreases slowly and shows residual order- 
ing up to large z distances. A comparison with previous 
DMC results on helium at the same density shows that 
the height of the ^He peaks is approximately an 80 % 
of the corresponding ones in II2, but with a z-ordering 
that also survives up to large distances; the additional 
20 % comes from the sizeable difference between the well 
depths of the respective interaction potentials. 

The significant differences in structure between the low 
and high density regimes are reflected even more clearly 



in the static structure factors. In Fig. 4, results for S{k) 
corresponding to the same densities reported in Fig. 3 
are shown. At high density, a characteristic result for a 
solid phase is obtained with a regular fc-spacing accord- 
ing to the only periodicity allowed by the ID geometry. 
At A = Ao, S{k) shows a first peak reflecting the local- 
ization observed in g(z) (Fig. 3), and a subsequent very 
smoothed maximum as expected in a homogeneous liquid 
phase. 

In conclusion, we have studied the zero-temperature 
equation of state of molecular hydrogen in a ID geometry 
and inside a narrow nanotube by means of the diffusion 
Monte Carlo method. The ID calculation predicts the 
existence of a self-bound system with a binding energy 
of -4.8 K and a quasi-continuous liquid-solid transition at 
high density. The comparison with a real system, hydro- 
gen in a nanotube, points to a close proximity between 
its properties and the ones of the ID limit. The predic- 
tion of a liquid II2 phase inside a (5,5) carbon nanotube 
is one of the main conclusions of the present work. The 
high one-dimensionality of this system would preclude a 
superfluid behavior but the use of wider nanotubes can 
provide a proper setup to observe the so-long desired su- 
perfluidity in molecular hydrogen. We expect the present 
work could encourage experimentalists to explore such a 
intriguing possibility. 
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FIG. 2. ID {Px, sohd hne) and T (P, dashed hne) pressures 
for H2 and '*He as a function of the linear density. 
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FIG. 1. Energy per particle of H2 as a function of the linear 
density. Open squares are the ID results, and filled squares 
are the T energies having subtracted the binding energy of a 
single molecule E),. The lines are the result of the polynomial 
fit (Eq. 4) with the optimal parameters reported in Table II. 
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FIG. 3. Two-body radial distribution function for ID H2 
at equilibrium (solid line) and at the liquid-solid transition 
density (dashed line). 
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FIG. 4. Static structure function for ID H2 near equilib- 
rium (open symbols) and at the liquid-solid transition density 
(filled symbols). 
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A (A~') E/N (ID, g = 0) E/N (ID, a ^ 0) E/N (T, a = 0) E/N (T, a ^ 0) 

0.329 98.083 ± 0.034 97.963 ± 0.016 -1453.99 ± 0.06 -1454.69 ± 0.04 

0.320 72.567 ± 0.013 72.523 ± 0.007 -1476.74 ± 0.05 -1476.88 ± 0.01 

0.312 53.264 ± 0.010 53.227 ± 0.010 -1493.790 ± 0.019 -1493.720 ± 0.002 

0.304 38.581 ± 0.018 38.636 ± 0.014 -1506.570 ± 0.03 -1506.540 ± 0.011 

0.290 19.203 ± 0.010 19.260 ± 0.003 -1523.730 ± 0.017 -1523.600 ± 0.02 



TABLE I. Energies per particle in K at high linear densities A for ID and T H2 systems, a = and a ^ correspond to the 
liquid and solid phases, respectively. 



Parameter ID H2 H2 in a tube 

Ao (A"') 0.2191 ± 0.0004 0.2200 ± 0.0006 

eo (K) -4.834 ± 0.007 -1544.880 ± 0.016 

A (K) 65.7 ± 3.6 69.7 ± 4.9 

B (K) 556.0 ± 46.6 429.7 ± 73.9 

X /v 1.98 1.5 

TABLE II. Parameters of the equation of state (Eq. 4) for the two systems studied. 
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